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Abstract 

Within the framework of QCD hght-cone sum rules (LCSR), we calculate the 
form factor for B^Dlv transitions with chiral current correlator. The resulting 
form factor depends on the distribution amplitude (DA) of the D-meson. We 
try to use three kinds of DA models of the D-meson. In the velocity transfer 
region 1.14 < y < 1.59, which renders the Operator Product Expansion (OPE) 
near light-cone x'^ = go effectively, the yielding behavior of form factor is in 
agreement with that extracted from the data on B ^ DW, within the error. In 
the large recoil region 1.35 < y < 1.59, the results are observed consistent with 
those of perturbative QCD (pQCD). The presented calculation can play a bridge 
role connecting those from the lattice QCD, heavy quark symmetry and pQCD to 
have an all-around understanding of B^Dlv transitions. 
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1 Introduction 



Calculation of the form factors for semileptonic transitions of B mesons has been being 
a subject discussed intensely. Recently, it has been shown that the B ^ tt transition 
form factor can be consistently analyzed by using the different approaches in the different 
q"^ regions [U El El HI • The perturbative QCD (pQCD) can be applied to the B ^ tt form 
factor in the large recoil (small g^) region and it is reliable when the involved energy scale 
is large enough [I]. The QCD light-cone sum rules (LCSR) can involve both the hard 
and soft contributions to the S — > tt form factor below ~ 18 GeV^ The lattice 
QCD simulations of the B ^ it transition form factor |3J are available only for the soft 
region q^ > 15 GeV^, because of the restriction to the vr energy smaller than the inverse 
lattice spacing. Thus the results from these three approaches might be complementary 
to each other. In Ref. |3] we recalculate the B n form factor in the pQCD approach, 
with the transverse momentum dependence included for both the hard scattering part 
and the nonperturbative wave functions(of vr and B) to get a more reliable pQCD result. 
By combining the results from these three methods we obtain a full understanding of the 
B ^ 71 transition form factor in its physical region < < {Mb — M^-)^ ^ 25 GeV^. 

Up to now, in comparison with heavy-to-light cases the calculations on heavy-to- 
heavy transitions can be done only for a certain specific kinematical range, although there 
have been a lot of discussions in the literature. In the BSW model j^I, the relevant form 
factors at zero momentum transfer are expressed as an overlap of initial and final meson 
wave functions for which they take the solutions of the Bethe-Salpeter (BS) equation in 
a relativistic harmonic oscillator potential. Then one extrapolates the result at = to 
the whole kinematical region assuming the nearest pole dominance. With the discovering 
of the heavy quark symmetry, the B^D form factor have been known better at zero 
recoil. This is because of the fact that in the heavy quark limit the resulting form 
factors — Isgur-Wise functions [H] at zero recoil are rigorously normalized. Including the 
leading symmetry breaking corrections, the deviation from this limit can be estimated 
at an order of a few percent due to Luke's theorem and therefore the value of the 
form factor at this point can be determined within a higher accuracy p . However, 
the dependence of the form factor on the velocity transfer y = v ■ v' (with v and v' 
being the velocities of the B and D mesons respectively) is difficult to get even in the 
leading order, in view of the arbitrary function (j{y) jH] which is introduced to simulate 
higher-resonances in the heavy quark effective theory (HQET). The lattice QCD, despite 
a rigorous nonperturbative approach, is just adequate to estimate the behavior of the 
form factors near the zero recoil 0. Among the other approaches are the QCD sum 
rules and pQCD. Ref.^ni applies the traditional 3-point sum rule to calculate the form 
factor at zero momentum transfer. It is concluded in Ref.[nJ that pQCD approach is 
applicable in the large recoil region and can give a consistent result with the experiment. 

It is necessary that there is a reliable estimate of B ^ D transition in the whole 
kinematically accessible range < g^ < (Mb—MdY ~ 11.6 GeV^, in order to account for 
the data on B^DW. For this purpose, it is practical, as shown in B ^ n case, to combine 
the result of QCD LCSR with those from the lattice QCD, heavy quark symmetry and 
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pQCD. The LCSR approach [12], where the non-perturbative dynamics are effectively 
parameterized in so-called light-cone wave functions, is regarded as an effective tool to 
deal with heavy-to-light exclusive processes. Although the B^D transition in question is 
a heavy-to-heavy one, the c-quark is much lighter compared to 6-quark and so discussing 
it with LCSR is plausible for the kinematical range where the OPE near light-cone = 
is valid. The other problem with our practical calculation is that the higher twist DA's 
of D meson, which are important but less studied, would enter into the sum rule result. 
However, an effective approach ^21 has been presented to avoid the pollution by some 
higher-twist DA's. This improved LSCR method uses a certain chiral current correlator 
as the starting point so that the relevant twist-3 wave functions make no contributions 
and the reliability of calculation can be enhanced to a large extend. Its applicability has 
been examined by a great deal of studies [21 El- In this paper we would like to employ 
the improved LCSR to discuss the form factor for the B—>-D transition and try to give 
a full understanding of QCD dynamics involved in the B — ^ Dlu. 

This paper is organized as follows. In the following Section we derive the LCSR for 
the form factor for 5 — D. A discussion of the DA models for the D-meson is given in 
section III. Section IV is devoted to the numerical analysis and comparison with other 
approaches. The last section is reserved for summary. 



2 Derivation of LCSR for The B^D Form Factor 

The B^D weak form factors f{q^) and f{q^) are usually defined as: 

<D{p)\c^l,h\B{p + q)> = 2f{q^)p^ + f{q^)q^, (1) 

with q being the momentum transfer. On the other hand, when applying the heavy 
quark symmetry to do discussion the following definition is advisable, 

<D{p)\c-i^h\B{p + q)> = ^niB mD[h+{y){v + v')^ + h_{y){v - v')^]. (2) 

If we neglect the masses of leptons in the decay final state of B—^DWi, only /(g^) is 
relevant and thus we can confine us to the discussion on f{q^)- Obviously, the following 
relation is observed between f{q^) and 

where J^B^oiy) = h+{y) - q'^ = ml + - 2mB niDy. 

Using the heavy quark symmetry, the value of form factor J^b^d{^) at zero recoil 
could be fixed better. Since in heavy quark limit = 1 and h-{l) = 0, the form 

factor JFb^d(I) should be close to 1. A systematic investigation gives JFb^d(I) = 
0.98 ± 0.07pn], with a less model dependence. This result is also confirmed with lattice 
calculations jH]. PQCD analyses are also made in the large recoil region y = 1.35 — 1.59, 
yielding a result consistent with the data. The LCSR calculation can help to understand 
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TB-*D{y) in the whole kinematical region in complementary to the lattice QCD with the 
heavy quark symmetry and pQCD approaches. 

To achieve a LCSR estimate of J-'B-^niy), we follow and use the following chiral 
current correlator n^(p, g): 



n^(p,g) = z J d''xe'P-<D{p)\T{cixh^{l + ^,)b{x),biO)t{l + ^,)d{0)}\0> 

= n(g2, (p + qf)p^ + Yi{q\ (p + qf)q^, (4) 

In the first place, we discuss the hadronic representation for the correlator. This can 
be done by inserting the complete intermediate states with the same quantum numbers 
as the current operator bi{l + 75)^. Isolating the pole contribution due to the lowest 
pseudoscalar B meson, we have the hadronic representation in the following: 

(p, q) = (p + qf)p^ + li'' {q\ {p + qf)q^ 

_ <D\c-f^b\B><B\bi-f5d\0> 
m% - {p + qY 

V rnl,-{p + qY ' 

Note that the intermediate states B^ contain not only the pseudoscalar resonance of 
masses greater than tub, but also the scalar resonances with = O"*", corresponding to 
the operator bd. With Eq.(^ and the definition of the decay constant J'b of the S-meson 

<B\bi-f5d\0> = mlfB/nib, (6) 

and expressing the contributions of higher resonances and continuum states in a form of 
dispersion integration, the invariant amplitudes 11^ and 11^ read, 

n^[g', (P + qf] = ^^(g')^^^^ + r _^^K^ds + subtractions, (7) 



mb(m| - {p + qY) J so s - {p + qY 
and 

n^'lg', (P + qf] = /^p'^^f^ + r —P^i^ds + subtractions, (8) 



mb{m% - {p + qY) J so s - {p + qY 

where the threshold parameter sq should be set near the squared mass of the lowest 
scalar B meson, the spectral densities p^{s) and p^{s) can be approximated by invoking 
the quark-hadron duality ansatz 

p''is){p''{s)) = p«^^(.)(p'3^^(.))^(. - .0). (9) 

On the other hand, we need to calculate the corrector in QCD theory to obtain the 
desired sum rule result. In fact, there is an effective kinematical region which makes OPE 
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applicable: (p + g)^ — mf-CO for the bd channel and q'^<{mb — rric)'^ — 2AQCD{^b — mc) 
for the momentum transfer. 

For the present purpose, it is sufficient to consider the invariant amplitude n(g^, (p + 
g)^) which contains the desired form factor. The leading contribution is derived easily 
by contracting the 6— quark operators to a free propagator: 

<o|r«|o>^/|^,e-l±^. (10) 

Substituting Eq.(fTm) into Eq.(jH), we have the two-particle contribution to the correlator, 

r (i^rd'^k 1 
nf ) = -2m,z J ^^e^(^-')^^^3^</^(p)|Tc(x)7M75rf(0)|0 >. (11) 

An important observation, as in Ref.^, is that only the leading non-local matrix ele- 
ment <i5(p)|c(x)7^75(i(0)|0> contributions to the correlator, while the nonlocal matrix 
elements <D{p)\c{x)i'y5d{0)\0> and <D(p)|c(x)cr^,y75(i(0)|0> whose leading terms are of 
twist 3, disappear from the sum rule. Proceeding to Eq.()ll|). we can expand the nonlocal 
matrix element <D(p)|Tc(a;)7^75(i(0)|0> as 

<D(p)|Tc(a;)7^75rf(0)|0> = -ip^^fo due'"'''"' ^ d{u) + higher twist terms, (12) 

where ipniu) is the twist-2 DA of D meson with u being the longitudinal momentum 
fraction carried by the c-quark, those DA's entering the higher-twist terms are of at least 
twist 4. The use of Eq. lfT^ yields 

n^*""?) \q^, (p + qf] = 2fDmb C du „ '^f^^^ — - + higher twist terms. (13) 

Jo mf — [up + qY 

Invoking a correction term due to the interaction of the b quark with a background 
field gluon into (10), the three-particle contribution 11^^'^^) is achievable. However, the 
practical calculation shows that the corresponding matrix element whose leading term is 
of twist 3 also vanishes. Thus, If we work to the twist-3 accuracy, only the leading twist 
DA is needed to yield a LCSR prediction. 

Furthermore, we carry out the subtraction procedure of the continuum spectrum, 
make the Borel transformations with respect to {p + q)^ in the hadronic and the QCD 
expressions, and then equate them. Finally, from Eq.(j2I) follows the LCSR for J^B^oiy), 
which is applicable to the velocity transfer region 1.14 < y < 1.59, 

2"^| jm^fD^ml/M^ 



/■I du 
X — exp 

J A U 



ml — (1 — M)(g^ — um'jj) 



^d{u), (14) 
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where 

{sq - - m|,)2 + Am\){ml - q^) - (sq - - "m]^) 
and = m'jj has been used 



A = ^^^^ — — (15) 



3 D-meson Distribution Amplitude 

Now let's do a brief discussion on an important nonperturbative parameter appearing 
in the LCSR formulalfT^. the leading twist DA of D-meson, ip£i{x). 

D-meson is composed of the heavy quark c and the light anti-quark q. The longitudi- 
nal momentum distribution should be asymmetry and the peak of the distribution should 
be approximately at x ~ 0.7. According to the definition in Eq. lfT^ . (poix) satisfies the 
normalization condition 



dxipD^x) = 1, (16) 

which is derived by the leptonic decay D — > fiu. 

In the pQCD calculations |llj, a simple model (we call model I) is adopted as 

ip^i\x) = Qxil - x){l - Cdil - 2x)) (17) 

which is based on the expansion of the Gegenbauer polynomials. Eq.()17p has a free 
parameter Cd which ranges from to 1. We will take = 0.7 as input. 

On the other hand, it was suggested in that the light-cone wave function of the 
Z^-meson be taken as: 



D 



+ 



X 



X 



(18) 



which is derived from the Brosky- Huang- Lepage (BHL) prescription jTTj. ^£)(x, k^) can 
be related to the DA by the definition 



(Pd{x) 



2v^ /-c/^k 
1^ 



(19) 



Substituting Eq. lfTH)) into Eq. lfT^ . we have a model of the DA(model II) 



D 



8772 52 



■x(l — x) exp 



D 



2 xm\ + (1 — x)m 



D 



xil — x) 



(20) 



where the parameters Ao and ho can be fixed by the normalization and the proba- 
bility of finding the \qq > Fock state in the D- meson, Pd 



D 



dx 



I^D(x,kx)|l 



(21) 
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As discussed in Ref . [TH] . Pd ~ 0.8 is a good approximation for the iI)-meson (As we have 
checked, change of Pd makes a numerical effect less than 2%). Then, taking Pd f» 0.8, 
fo = 240MeV, = 1.3GeV and = 0.35GeV, we have Ad = 63.6GeV-\ bl = 
0.292GeV-l 

Furthermore, as argued in Ref.jTHj, a more complete form of the light-cone wave 
function should include the Melosh rotation effect in spin space: 



x,kj 



XD{x,k±)exp 



-hi 



'k^ + ml 



X 



kij-m|^ 
1 — X 



(22) 



with the Melosh factor, 



Xd{x, kx) 



(1 — x)mc + xrud 
k^ + ((1 — x)mc + xmd) 



(23) 



It can be seen from Eq. (j23p that XD{Xik_[_) — + 1 as tjIc —>■ oo, since there is no spin 
interaction between the two quarks in the heavy-flavor meson, ie., the spin of the heavy 
constituent decouples from the gluon field, in the heavy quark limit jH]. However the 
c-quark is not heavy enough to neglect the Melosh factor. 

After integration over k^ the full form of D meson DA can be achieved (model III): 



Ad^3x{1-x) 



I -Erf 



boy 



'x(l — x) 



exp 



(xm^ + {I — x)m^ — y"^) 



x(l — x) 



(24) 

where y = xrrid + (1 — x)mc and the error function Erf{x) is defined as Erf{x) = 
- Jq exp{—t'^)dt. Using the same constraints as in Eq.lfTBj) and (EH), the parameters Ajy 
and bo are fixed as Ad = 62.8GeV-^ and b^ = 0.265GeV~^ 

In this paper we will employ the above three kinds of models to do numerical cal- 
culation. All these DA's of the Z)-meson are plotted in Fig.([[]) for a comparison. It is 
shown that they are of similar shape and all of them exhibit a maximum at x ~ 0.6 — 0.7 
as expected. 



4 Numerical Result and Discussion 

Apart from the DA of D-meson, the decay constant of 5-meson fs is among the 
important nonperturbative inputs. For consistency, we use the following corrector 

K{q^)=t I d^xe''^^ < 0|g(x)(l +75)&(a^),&(0)(l -75)g(0)|0 >, (25) 

to recalculate it in the two-point sum rules. The calculation should be limited to leading 
order in QCD, since the QCD radiative corrections to the sum rule for J-'B^oiy) are not 
taken into account. The value of the threshold parameter s'q is determined by a best fit 
requirement in the region 10GeV^<M^<20GeV^, where is the corresponding Borel 
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nib 


s'o Ib 


So 


set 1 
set 2 
set 3 


4.85 
4.80 
4.75 


29.5 0.076 
29.8 0.090 
30.0 0.103 


30.3 
30.8 
31.3 



Table 1: Parameter sets for fs and J-'B-^niy), Sq and sq for Jb and Ts^oiy) respectively; 
nib and /b are given in GeV, sq and s'q in GeV^. 

parameter. The same procedure is performed for J^B^niy), resulting in different values 
of the threshold parameter sq. The result is listed in Tab^ Choosing sq and s'q in this 
way, the dependence of fs and J-'B^oiy) on the Borel parameter is very weak and thus 
we can simply evaluate them at = = 15GeV^. The other input parameters 
are taken as = 5.279GeV, m^) = 1.869GeV. As we have ignored all the radiation 
corrections, we don't expect our values of to be good predictions of that quantity. 

With the parameters chosen, it is straightforward to calculate the form factor J^s^oiy) 
in the region 1.14 < y < 1.59. The results with different sets of parameters are plotted 
in Fig.Q, where only model II has been used for simplicity. It is shown that the change 
of parameters can induce a uncertainty of about 10 — 15% if we let vary between 
4.75 — 4.85 GeV. By fitting the data, the behavior of Ts^Diy) has been known using 
the parametrization 

rB^D{y) = ^b^d(1)[1 - pUy - 1) + cniy - \f + 0((y - 1)^)], (26) 

with 

= 0.69 ±0.14, cz) = 0, 

= 0.69^°:^^, = 0.00^°:^^, (27) 

corresponding to the linear and quadratic fits ^H], respectively. With the three DA 
models, the resulting dependence of TB^oiy) on the velocity transfer y, along with 
that extracted experimentally is illustrated in Fig.Q. In what follows, we denote the 
LCSR results for the form factor by S'lid those extracted experimentally by 

T^b^dW)- For comparison, a figure-copy which expresses the pQCD results in |TTj is 
given in Fig.Q. In the region to which the LCSR method is applicable, the central 
values of J^B^j^iy) turn out to be a bit smaller than the corresponding those of J^^ZDiy)^ 
using the DA models II and III as inputs; however, both of them are in accordance 
with each other within the error. The situation with model I DA is about the same. 
The central value of the form factor at the largest recoil is ^-"^^£,(1.59) = 0.58 versus 
^b££)(1-59) = 0.40 — 0.50, depending on the DA models. We note that the behavior 
of J^B^^iy) is essentially unchanged when the three different DA's are used. From the 
present calculations, therefore it is too early to draw a conclusion which DA model is 
more suitable to reflect the characteristics of QCD dynamics inside the D meson. When 
a comparison is made between the pQCD and LCSR predictions, the consistent results 
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can also be observed at the larger recoil. Of course, the two approaches describe the 
different dynamics in B ^ D transitions. Whereas the use of LCSR approach is to 
assume that the soft exchanges dominate in the weak decay in question, applying pQCD 
method to do calculation corresponds to the viewpoint that the hard exchanges do. 

In fact, the kinematical region we give, which makes LCSR results valid, is a con- 
servative estimate. It is possible to extrapolate the present LCSR calculation to the 
small recoil region. If it is true, we find that in the whole kinematically accessible 
range 1.0 < y < 1.59, the yielding LCSR estimates are compatible with the data. 
For instance, at zero recoil it follows that — 1-02 (using model III), which 

is in a good agreement with the evaluation obtained using the heavy quark symmetry: 
J-b^d{^) = 0.98 ± 0.07. Nevertheless, we have to emphasize that a full understanding 
of the dynamics involved in B ^ D transition should be obtained by combining the 
three different approaches — the lattice QCD calculations with the heavy quark symme- 
try considered, LCSR results and pQCD predictions, which are complementary to each 
other. The LCSR results with chiral current correlator may act as a bridge connecting 
those of other approaches. 

5 Summary 

We have discussed the form factor for B^D transitions TB^oiy)-, using the improved 
QCD LCSR approach where with the chiral current correlator chosen only the leading 
twist DA of the Z^-meson is relevant at twist-3 accuracy. The resulting LCSR's for 
^B^oiy) are available in the velocity transfer region 1.14 < y < 1.59. Calculation is 
done using three different twist-2 DA models for D meson. It is shown the numerical 
results are less sensitive to the choice of DA, and are of a central value slight smaller 
than but within the error in a agreement with those obtained by fitting the data on 
B —>■ DW. In the larger recoil region 1.35 < y < 1.59 where pQCD is applicable, the 
results presented here are consistent with ones of pQCD. From the practical calculations, 
we find that the present results might be extrapolated to the smaller recoil region so that 
the B^D transitions are calculable in the whole kinematically accessible range, using 
the improved LCSR approach. 

Also, we argue that for understanding the form factor for B^DW in the whole kine- 
matical range a combined use is necessary of three different methods: the lattice QCD 
(with the heavy quark symmetry considered), improved LCSR and pQCD approaches, 
which are adequate to do calculation in different kinematical regions and so could be 
complementary to each other. The LCSR approach plays a bridge role in doing such 
calculation. 

The present findings can be improved once the QCD radiative correction to the LCSR 
is taken into account and a more reliable twist-2 DA of D meson becomes available. From 
the previous discussion in [12] , however, it is expected that the QCD radiative correction 
can not change the present results too much. 
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Figure 1: Different kinds of D-meson DAs, solid and dashed curves correspond to model 
III and II, while the dotted hne expresses model I. 




Figure 2: Dependence of Tb^d on the different sets of parameters rub, fB,so. The three 
curves correspond to the the parameter set 1—3 from bottom to top. Here we use model 
II for the D-meson DA for simplicity. 
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Figure 3: Tb^d as a function of the velocity transfer (with the parameters in the set 2). 
The thin hues expresses the experiment fits results, the solid line represents the central 
values, the dashed(dash-dotted) lines give the bounds from the linear(quadratic) fits. 
The thick lines correspond to our results, with the solid, dashed and dash-dotted lines 
for model III, II and I respectively. 



form factor 
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(a) 



Figure 4: pQCD results for J^B^niy) copied from Ref . [TT] . As in Fig.(jni), the solid line 
represents the central values, the dashed(dash-dotted) hues give the bounds from the 
linear(quadratic) fits. The circles corresponds to pQCD results using model I for the 
D-meson DA with Co = 0.5, 0.7, 0.9 from bottom to top. 
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